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The 7-ray burst (GRB) model for production of ultra-high-energy, > 10 19 eV, cosmic-rays is based on the 
hypothesis that GRBs arise from the dissipation of the kinetic energy of relativistic fireballs at cosmological 
distances. Recent observations of delayed low energy emission, "afterglow," from GRB sources strongly support 
the validity of this hypothesis. Observations also provide quantitative support for the model. The inferred physical 
fireball parameters imply that protons may be accelerated to > 10 20 eV, and the inferred GRB energy generation 
rate is similar to that required to account for the observed flux of ultra- high-energy cosmic-rays (UHECRs). 

Strong suppression of cosmic-ray flux is expected in this model above 10 19 ' 7 eV, due to proton interaction with 
microwave background photons. Strong deviations from model flux derived under the assumption of uniform 
source distribution is expected above 10 20 eV, due to source discreteness and due to inhomogeneities in source 
distribution. In particular, the flux above 10 20 ' 5 eV is expected to be dominated by few, narrow spectrum sources. 
While model predictions can not be tested (with high confidence level) using present data, the predicted signatures 
should be observed with the planned Auger and Telescope-Array UHECR detectors. 

A natural consequence of the GRB model of UHECR production is the conversion of a large fraction, ~ 10%, 
of the fireball energy to accompanying burst of ~ 10 14 eV and ~ 10 18 eV neutrinos. A km 2 neutrino detector 
would observe several tens of events per year correlated with GRBs, and test for neutrino properties (e.g. flavor 
oscillations, for which upward moving t's would be a unique signature, and coupling to gravity) with an accuracy 
many orders of magnitude better than is currently possible. 



1. Introduction 

The origin of GRBs, bursts of 0. 1 MeV— 1 MeV 
photons lasting for a few seconds, remained un- 
known for over 20 years, primarily because GRBs 
were not detected prior to 1997 at wave-bands 
other than 7-rays (see pj| for review of 7-ray ob- 
servations). The isotropic distribution of bursts 
over the sky suggested that GRB sources lie at 
cosmological distances, and general phenomeno- 
logical considerations were used to argue that the 
bursts are produced by the dissipation of the ki- 
netic energy of a relativistic expanding fireball 
(see H for review). 

Adopting the cosmological fireball hypothesis, 
it was shown that the physical conditions in the 
fireball dissipation region allow Fermi accelera- 
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tion of protons to energy > 10 20 eV ||@], and that 
the average rate at which energy is emitted as 7- 
rays by GRBs is comparable to the energy gen- 
eration rate of UHECRs in a model where UHE- 
CRs are produced by a cosmological distribution 
of sources ||. Based on these two facts, it was 
suggested that GRBs and UHECRs have a com- 
mon origin (see || for review). 

In the last two years, afterglows of GRBs have 
been discovered in X-ray, optical, and radio wave 
bands (see JtJ for review). Afterglow observa- 
tions confirmed the cosmological origin of the 
bursts, through the rcdshift determination of sev- 
eral GRB host-galaxies (see || for an updated 
list), and confirmed || standard model predic- 
tions jl(| of afterglows that result from the colli- 
sion of an expanding fireball with its surrounding 
medium. These observations therefore provide 
strong support for the GRB model of UHECR 
production. 

In this review, UHECR and neutrino produc- 
tion in GRBs is discussed in the light of recent 
GRB and UHECR observations. The fireball 
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model is briefly described in §2.1, and proton ac- 
celeration in GRB fireballs is discussed in §2.2. 
Recent claims, according to which protons can 
not be accelerated to > 10 20 eV in the fireball 
fill , are shown in §2.2 to be erroneous. Implica- 
tions of recent afterglow observations to high en- 
ergy particle production are discussed in §3. It is 
shown that, contrary to some recent claims |12| , 
the GRB energy generation rate implied by af- 
terglow observations is similar to the energy gen- 
eration rate required to account for the flux of 
> 10 19 eV cosmic-rays. Model predictions are 
shown to be consistent with the observed UHECR 
spectrum in §4. 

Predictions of the GRB model for UHECR pro- 
duction, that can be tested with future UHECR 
experiments, are discussed in §5. Implications of 
the detection by the AGASA experiment of mul- 
tiple high energy events with consistent arrival di- 
rections 13 1 is also discussed in §5. High energy 



neutrino production in fireballs and its implica- 
tions for future high energy neutrino detectors are 
discussed in §6. 

2. UHECR from GRB fireballs 

2.1. The fireball model 

In the fireball model of GRBs (m) , a compact 
source, of linear scale r ~ 10 7 cm, produces 
a wind characterized by an average luminosity 
L ~ 10 52 ergs _1 and mass loss rate M = L/rjc 2 . 
At small radius, the wind bulk Lorentz factor, T, 
grows linearly with radius, until most of the wind 
energy is converted to kinetic energy and T sat- 
urates at T ~ i] ~ 300. Variability of the source 
on time scale At, resulting in fluctuations in the 
wind bulk Lorentz factor T on similar time scale, 
then leads to internal shocks jlj] in the expanding 
fireball at a radius 



T 2 cAt = 3 x 10 1 



T300^10ms cm, 



(1) 



where F = 300r 300 , At = 10Aii Oms ms. If the 
Lorentz factor variability within the wind is sig- 
nificant, internal shocks would reconvert a sub- 
stantial part of the kinetic energy to internal 
energy. It is assumed that this energy is then 
radiated as 7-rays by synchrotron and inverse- 
Compton emission of shock-accelerated electrons. 



In this model, the observed 7-ray variability 
time, ~ ri/T 2 c ps At, reflects the variability time 
of the underlying source, and the GRB duration, 
T ~ 10s, reflects the duration over which energy 
is emitted from the source. A large fraction of 
bursts detected by BATSE show variability on 
the shortest resolved time scale, ~ 10 ms JlfJ, 
and some show variability on shorter time scales, 
~ 1 ms [0. This sets the constraint on under- 
lying source size, tq < cAt ~ 10 7 cm. The wind 
must be expanding relativistically, with a Lorentz 
factor r ~ 300, in order that the fireball pair- 
production optical depth be small for observed 
high energy, ~ 100 MeV, GRB photons @. 

The wind Lorentz factor is expected to fluctu- 
ate on time scales ranging from the source dy- 
namical time, At, to the wind duration T, lead- 
ing to internal collisions over a range of radii, 



n = T 2 cAt to 



T 2 cT. Internal shocks 



are generally expected to be "mildly" relativistic 
in the fireball rest frame, i.e. characterized by 
Lorentz factor 74 — 1 ~ 1, since adjacent shells 
within the wind are expected to expand with 
Lorentz factors which do not differ by more than 
an order of magnitude. 

As the fireball expands, it drives a relativistic 
shock (blastwave) into the surrounding gas. At 
early time, the fireball is little affected by this 
external interaction. At late time, most of the 
fireball energy is transferred to the surrounding 
gas, and the flow approaches self-similar expan- 
sion. For typical fireball parameters, the transi- 
tion to self-similar expansion occurs at a radius 
r oj T 2 cT. At this radius, mildly relativistic re- 
verse shocks propagate into the fireball ejecta and 
decelerate it Jl{|. The reverse shocks disappear 
on (observed) times scale T, and the flow becomes 
self-similar at later time, with a single, relativis- 
tic decelerating shock propagating into the sur- 
rounding medium. Plasma conditions in the re- 
verse shocks are similar to those of internal shocks 
arising from variability on time scale ~ T, since 
both are mildly relativistic and occur at similar 
radii. In the discussion that follows we therefore 
do not discuss the reverse shocks separately from 
the internal shocks. 

The shock driven into the ambient medium 
continuously heats new gas, and accelerates rel- 
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ativistic electrons that may produce by syn- 
chrotron emission the delayed radiation, "after- 
glow," observed on time scales of days to months. 
As the shock-wave decelerates, the emission shifts 
with time to lower frequency. 

2.2. Fermi acceleration in GRBs 

In the fireball model, the observed GRB and 
afterglow radiation is produced by synchrotron 
emission of shock accelerated electrons. In the 
region where electrons are accelerated, protons 
are also expected to be shock accelerated. This 
is similar to what is thought to occur in super- 
novae remnant shocks |p0(| . We consider below 
proton acceleration in internal (and reverse) fire- 
ball shocks. Since the internal shocks are mildly 
relativistic, we expect results related to particle 
acceleration in sub-relativistic shocks (see, e.g., 
20 1 for review) to be valid for the present sce- 
nario. In particular, the predicted energy distri- 
bution of accelerated protons is dN p /dE p oc E~ 2 . 

Two constraints must be satisfied by fireball 
wind parameters in order to allow proton accel- 
eration to E p > 10 20 eV in internal shocks S: 

W£e > °- 02r 300 £; p,20 i 7, 1 52' 



(2) 



and 

r > lSO^o^Aiio 1 ^. 



(3) 



Here, E p = 10 20 E p ^o eV, L 7 = 10 52 L 7j52 erg/s 
is the 7-ray luminosity, £s is the fraction of the 
wind energy density which is carried by magnetic 
field, 47rr 2 cr 2 (S 2 /87r) = £ B L, and £ e is the frac- 
tion of wind energy carried by shock accelerated 
electrons. Since the electron synchrotron cool- 
ing time is short compared to the wind expansion 
time, electrons lose their energy radiatively and 
L is related to the observed 7-ray luminosity by 
L 7 « £ e L. The first condition must be satisfied 
in order for the proton acceleration time t a to be 
smaller than the wind expansion time. The sec- 
ond condition must be satisfied in order for the 
synchrotron energy loss time of the proton to be 
larger than t a . 

From Eqs. (||) and (||), we infer that a dissipa- 
tive ultra-relativistic wind, with luminosity and 
variability time implied by GRB observations, 



satisfies the constraints necessary to allow the ac- 
celeration of protons to energy > 10 20 eV, pro- 
vided that the wind bulk Lorentz factor is large 
enough, T > 100, and that the magnetic field 
is close to equipartition with electrons. The for- 
mer condition, T > 100, is remarkably similar to 
that inferred based on 7-ray spectra. There is 
no theory at present that allows a basic princi- 
ples calculation of the strength of the magnetic 
field. However, magnetic field close to equiparti- 
tion, £g ~ 1, is required in order to account for 
the observed 7-ray emission (see also §3). 

We have assumed in the discussion so far that 
the fireball is spherically symmetric. However, 
since a jet-like fireball behaves as if it were a 
conical section of a spherical fireball as long as 
the jet opening angle is larger than T -1 , our re- 
sults apply also for a jet-like fireball (we are inter- 
ested only in processes that occur when the wind 
is ultra-relativistic, T ~ 300, prior to significant 
fireball deceleration). For a jet-like wind, L in 
our equations should be understood as the lumi- 
nosity the fireball would have carried had it been 
spherically symmetric. 

It has recently been pointed out in [[n] that 
conditions at the external, highly relativistic 
shock driven by the fireball into the ambient 
gas are not likely to allow proton acceleration to 
ultra- high energy. Although correct, this observa- 
tion is irrelevant to the scenario considered here 
based on ||, since in this scenario protons are 
accelerated in internal, mildly relativistic fireball 
shocks. 

3. Implications of afterglow observations 

In addition to providing support to the valid- 
ity of the qualitative fireball scenario described in 
§2.1 0, afterglow observations provide quantita- 
tive constraints on fireball model parameters. 

The determination of GRB redshifts implies 
that the characteristic GRB 7-ray luminosity and 
emitted energy, in the 0.05 to 2 MeV band, 
are L 1 ~ 10 52 erg/s and E 1 ~ 10 53 erg respec- 
tively (e.g. J8|), an order of magnitude higher 
than the values assumed prior to afterglow de- 
tection (here, and throughout the paper, we as- 
sume an open universe, f2 = 0.2, A = 0, and 
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Hq = 75 km/s Mpc). The increased GRB lumi- 
nosity scale implies that the constraint @ on the 
fireball magnetic field is less stringent than pre- 
viously assumed. 

Due to present technical limitations of the ex- 
periments, afterglow radiation is observed in most 
cases only on time scale >> 10 s. At this stage, 
radiation is produced by the external shock driven 
into the surrounding gas, and afterglow observa- 
tions therefore do not provide direct constraints 
on the magnetic field energy fraction £g at the 
internal and reverse shocks, where protons are ac- 
celerated to ultra-high energy. In one case, how- 
ever, that of GRB 990123, reverse shock emission 
was detected over ~ 10 s time scal e [[21||22| ]. For 
this case, the inferred value of £b J23[ is consis- 
tent with the constraint (Q). Clearly, more ob- 
servations are required to determine whether this 
condition is generally satisfied. 

The observed GRB redshift distribution im- 
plies a GRB rate of -Rgrb ~ 10/Gpc 3 yr at z <~ 1. 
The present, z = 0, rate is less well constrained, 
since most observed GRBs originate at redshifts 
1 < z < 2.5 Present data are consis- 

tent with both no evolution of GRB rate with 
redshift, and with strong evolution (following, 
e.g., the luminosity density evolution of QSOs or 
the evolution of star formation rate), in which 

r G rb(z = i)/Rgbb(z = °) ~ 8 & The ener sy 

observed in 7-rays reflect the fireball energy in 
accelerated electrons. Afterglow observations im- 
ply that accelerated electrons and protons carry 
similar energy ||. Thus, the inferred z = rate 
of cosmic-ray production by GRBs is similar to 
the generation rate of 7-ray energy, 

E 2 {dh CR /dE) z=a = 10 44 Cerg/Mpc 3 yr, (4) 

where £ is in the range of ~ 1 to ~ 8. This 
energy generation rate is remarkably similar to 
that implied by the observed UHECR flux (see 

4 It has recently been argued Jl2| that the 2 = 
GRB 7-ray energy generation rate is much smaller, ~ 
10 42 erg/Mpc 3 yr. Most of the discrepancy between this 
result and our result can be accounted—for by noting two 
errors made in the analysis of ref. ]12j ]: estimating the 
energy generation rate as the product of the GRB rate 
and the median, rather than average, GRB energy, and 




Figure 1. The UHECR flux expected in a cos- 
mological model, where high-energy protons are 
produced at a rate (E 2 dhcR/dE) z=0 = 0.8 x 
10 44 erg/Mpc 3 yr as predicted in the GRB model 
[Eq. (||)], compared to the Fly's Eye, Yakutsk 
and AGASA data, la flux error bars are shown. 
The highest energy points are derived assuming 
the detected events (1 for Fly's Eye and Yakutsk, 
4 for AGASA) represent a uniform flux over the 
energy range 10 20 eV-3 x 10 20 eV. 



4. Comparison with UHECR observations 

Fly's Eye Q and AGASA g@ results con- 
firm the flattening of the cosmic-ray spectrum at 
~ 10 19 eV, evidence for which existed in previous 
experiments with weaker statistics P9|| . Fly's Eye 
data is well fitted in the energy range 10 17 6 eV to 
10 19 ' 6 eV by a sum of two power laws: A steeper 
component, with differential number spectrum 
J oc E~ 3 - 50 , dominating at lower energy, and a 
shallower component, J oc E~ 2 ' 61 , dominating at 
higher energy, E > 10 19 eV. The flattening of the 
spectrum, combined with the lack of anisotropy 
and the evidence for a change in composition from 
heavy nuclei at low energy to light nuclei (pro- 
tons) at high energy |3(J, suggest that an extra- 

using (following EHj) the GRB energy observed in the 50 
to 300 keV band, where only a small fraction of the 0.05 
to 2 MeV 7-ray energy is observed. 
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Galactic source of protons dominates the flux at 
high energy E > 1 19 eV. 

In Fig. 1 we compare the UHECR spectrum, 
reported by the Fly's Eye @, the Yakutsk [|lj, 
and the AG AS A |28| experiments, with that pre- 
dicted by the GRB model. The proton genera- 
tion rate is assumed to evolve in redshift follow- 
ing QSO luminosity evolution |}2| . Note, that the 
cosmic-ray spectrum at energy > 10 19 cV is little 
affected by modifications of the cosmological pa- 
rameters or of the redshift evolution of cosmic-ray 
generation rate, since cosmic-rays at this energy 
originate from distances shorter than several hun- 
dred Mpc. The spectrum and flux at E > 10 19 eV 
is mainly determined by the present (z — 0) gen- 
eration rate and spectrum. The absolute flux 
measured at 3 x 10 18 eV differs between the var- 
ious experiments, corresponding to a systematic 
~ 10% (~ 20%) over-estimate of event energies in 
the AGASA (Yakutsk) experiment compared to 
the Fly's Eye experiment (see also [0). In Fig. 
1, the Yakutsk energy normalization is used. 

The suppression of model flux above 10 19 7 eV 
is due to energy loss of high energy protons in 
interaction with the microwave background, i.e. 
to the "GZK cutoff" f|. Both Fly's Eye and 
Yakutsk data show a deficit in the number of 
events, compared to the number expected based 
on extrapolation of the J oc E~~ 2 61 power-law fit, 
consistent with the predicted suppression. The 
deficit is, however, only at a 2er confidence level 
|U . The AGASA data is consistent with Fly's Eye 
and Yakutsk results below 10 20 eV. A discrepancy 
may be emerging at higher energy, > 10 20 eV, 
where the Fly's Eye and Yakutsk experiments de- 
tect 1 event each, and the AGASA experiment 
detects 6 events for similar exposure. 

The flux above 10 20 eV is dominated by sources 
at distances < 30 Mpc || (see §5). Since 
the distribution of known astrophysical systems 
(e.g. galaxies, clusters of galaxies) is inhomoge- 
neous on scales of tens of Mpc, significant devi- 
ations from model predictions presented in Fig. 
1 for a uniform source distribution are expected 
above 10 20 eV. It has recently been shown [^4| 
that clustering of cosmic-ray sources leads to a 
standard deviation, er, in the expected number, 
N, of events above 10 20 eV, given by a/N = 



0.9(<io/10Mpc) ' 9 , where do is the unknown scale 
length of the source correlation function and cfo ~ 
10 Mpc for field galaxies. 

An order of magnitude increase in the exposure 
of UHECR experiments, compared to that avail- 
able at present, is required to test for the exis- 
tence of the GZK cutoff || . Such exposure would 
allow this test through an accurate determination 
of the spectrum in the energy range of 10 19 7 eV 
to 10 20 eV, where the effects of source inhomo- 
geneities are expected to be small H|||. More- 
over, an order of magnitude increase in exposure 
will also allow to determine the source correlation 
length do, through the detection of anisotropies in 
the arrival directions of ~ 10 19 5 eV cosmic-rays 
over angular scales of 9 ~ <io/30 Mpc ]34| . 

Finally, we note that preliminary results from 
the HiRes experiment were presented in this con- 
ference |35|, reporting 7 events beyond 10 20 eV for 
an exposure similar to that of the Fly's Eye. It 
is difficult to decide how to interpret this result, 
since the discrepancy between HiRes and Fly's 
Eye results is present not only above 10 20 eV but 
also at lower energy, where Fly's Eye, AGASA 
and Yakutsk experiments are in agreement: 13 
events above 6 x 10 19 eV are reported in the pre- 
liminary HiRes analysis, while only 5 events at 
that energy range are reported by the Fly's Eye. 
We therefore believe that unambiguous conclu- 
sions based on the recent HiRes data can only be 
drawn after a complete analysis of the HiRes data 
is published. 

5. GRB model predictions for planned 
UHECR experiments 

The energy of the most energetic cosmic ray 
detected by the Fly's Eye experiment is in ex- 
cess of 2 x 10 20 eV, and that of the most energetic 
AGASA event is ~ 2 x 10 20 eV. On a cosmo- 
logical scale, the distance traveled by such ener- 
getic particles is small: < lOOMpc (50Mpc) for 
the AGASA (Fly's Eye) event (e.g., |§). Thus, 
the detection of these events over a ~ 5yr period 
can be reconciled with the rate of nearby GRBs, 
~ 1 per 100 yr out to lOOMpc, only if there is a 
large dispersion, > lOOyr, in the arrival time of 
protons produced in a single burst. 
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The required dispersion is likely to occur due 
to the combined effects of deflection by random 
magnetic fields and energy dispersion of the par- 
ticles |U . A proton of energy E propagating over 
a distance D through a magnetic field of strength 
B and correlation length A is deflected by an an- 
gle 9 S ~ (D/ X) 1 ^ 2 X/ Rl, which results in a time 
delay, compared to propagation along a straight 
line, t(E,D) ps 9 2 s D/Ac oc B 2 X. The random en- 
ergy loss UHECRs suffer as they propagate, ow- 
ing to the production of pions, implies that at 
any distance from the observer there is some fi- 
nite spread in the energies of UHECRs that are 
observed with a given fixed energy. For protons 
with energies > 10 20 eV the fractional RMS en- 
ergy spread is of order unity over propagation 
distances in the range 10 — lOOMpc (e.g. p6[). 
Since the time delay is sensitive to the particle 
energy, this implies that the spread in arrival time 
of UHECRs with given observed energy is com- 
parable to the average time delay at that energy, 
t(E,D). 

The magnetic field required in order to pro- 
duce a spread r(E = 10 20 eV, D = lOOMpc) > 
100 yr, is well below the current upper bound 
on the inter-galactic magnetic field, BX 1 ^ 2 < 
10~ 9 G Mpc 1/2 @, which allows a spread r ~ 
10 5 yr. We discuss below some implications, 
unique to the GRB model, of time delays induced 
by magnetic-field deflection. 

5.1. The highest energy sources 

The rapid increase with energy of the pion 
production energy loss rate effectively introduces 
a cutoff distance, D C (E), beyond which sources 
do not contribute to the flux above E. The 
function D C (E) is shown in Fig. 2. We de- 
fine a critical energy E c , for which the aver- 
age number of sources at D < D c (Eg) is f, 



-RgrbD c (E c ) 3 t[E c ,D c (E c )] = 1 



Al- 



though E c depends through r on the unknown 
properties of the intergalactic magnetic field, the 
rapid decrease of D C (E) with energy near 10 20 eV 
implies that E c is only weakly dependent on the 
value of B 2 X. In The GRB model, the prod- 
uct R grb t(D = lOOMpc, £ = 10 20 eV) is ap- 
proximately limited to the range 10~ 6 Mpc 
to 10~ 3 Mpc~ 3 (The lower limit is set by the 
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Figure 2. Results of a Monte-Carlo realization 
of the bursting sources model, with E c = 1.4 x 
10 20 eV: Thick solid line- overall spectrum in 
the realization; Thin solid line- average spectrum, 
this curve also gives D C (E); Dotted lines- spectra 
of brightest sources at different energies. 



requirement that at least a few GRB sources 
be present at D < 100 Mpc, and the upper 
limit by the Faraday rotation bound BX 1 / 2 < 
10~ 9 G Mpc 1/2 and R GRB < 10/ Gpc 3 yr)- 
The corresponding range of values of E c is 
10 20 eV < E c < 3 x 10 20 eV. 

Fig. 2 presents the flux obtained in one re- 
alization of a Monte-Carlo simulation described 
in ref. (3£| of the total number of UHECRs re- 
ceived from GRBs at some fixed time for E c = 
1.4 x 10 20 eV. For each realization the distances 
and times at which cosmological GRBs occurred 
were randomly drawn. Most of the realizations 
gave an overall spectrum similar to that presented 
in Fig. H when the brightest source of this real- 
ization (dominating at 10 20 eV) is not included. 
At E < E c , the number of sources contributing 
to the flux is very large, and the overall UHECR 
flux received at any given time is near the average 
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flux (obtained for spatially and temporally homo- 
geneous UHECR volume emissivity). At E > E c , 
the flux will generally be much lower than the 
average, because there will be no burst within a 
distance D C (E) having taken place sufficiently re- 
cently. There is, however, a significant probabil- 
ity to observe one source with a flux higher than 
the average. A source similar to the brightest one 
in Fig. H appears ~ 5% of the time. 

At any fixed time a given burst is observed 
in UHECRs only over a narrow range of energy, 
because if a burst is currently observed at some 
energy E then UHECRs of much lower (higher) 
energy from this burst will arrive (have arrived) 
mainly in the future (past). For energies above 
the pion production threshold, E > 10 19 - 7 eV, the 
dispersion in arrival times of UHECRs with fixed 
observed energy is comparable to the average de- 
lay at that energy. This implies that the spec- 
tral width AE of the source at a given time is 
of order the average observed energy, AE ~ E. 
Thus, bursting UHECR sources should have nar- 
rowly peaked energy spectra, and the brightest 
sources should be different at different energies. 
For steady state sources, on the other hand, the 
brightest source at high energies should also be 
the brightest one at low energies, its fractional 
contribution to the overall flux decreasing to low 
energy only as D c (E)~ l . A detailed numerical 
analysis of the time dependent energy spectrum 
of bursting sources is given in |^| . 

The AGASA experiment reported the presence 
of one triplet and three doublets of UHECRs 
with angular separations (within each multiplet) 
< 2.5°, roughly consistent with the measure- 
ment error, among a total of 47 UHECRs with 
E > 4 x 10 19 eV U). The probability to have 
found such multiplets by chance is ~ 1%. There- 
fore, this observation favors the bursting source 
model, although more data are needed to confirm 
it. 

Testing the GRB model predictions described 
above requires an exposure 10 times larger than 
that of present experiments. Such increase is ex- 
pected to be provided by the planned Auger jioj 
and Telescope Array |44[] detectors. 



5.2. Spectra of Sources at E < 4 x 10 19 eV 

For nearby, D < 100 Mpc, sources contributing 
at E < 4 x 10 19 eV, pion production energy loss is 
negligible, and particle energy may be considered 
constant along the propagation path. In this case, 
the spectral shape of individual sources depends 
primarily on the magnetic field correlation length 

If A > D6 S (D,E) ~ D(D/Xy/ 2 X/R L , all 
UHECRs that arrive at the observer are essen- 
tially deflected by the same magnetic field struc- 
tures, and the absence of random energy loss 
during propagation implies that all rays with a 
fixed observed energy would reach the observer 
with exactly the same direction and time de- 
lay. At a fixed time, therefore, the source would 
appear mono-energetic and point-like (In real- 
ity, energy loss due to pair production results 
in a finite but small spectral and angular width, 

ae/e ~ se/e s < i% @). 

If, on the other hand, A < D9 S (D, E), the de- 
flection of different UHECRs arriving at the ob- 
server are essentially independent. Even in the 
absence of any energy loss there are many paths 
from the source to the observer for UHECRs of 
fixed energy E that are emitted from the source 
at an angle < 8 S relative to the source-observer 
line of sight. Along each of the paths, UHECRs 
are deflected by independent magnetic field struc- 
tures. Thus, the source angular size would be of 
order 9 S and the spread in arrival times would be 
comparable to the characteristic delay r, leading 
to AE/E ~ 1 (The spectral shape of sources is 
given in analytic form for this case in J42fl). 

For D = 30Mpc and E ~ 10 19 eV, the 8 S D = X 
line divides the allowed region (for the GRB 
model) in the B-X plane at A ~ IMpc. Thus, 
measuring the spectral width of bright sources 
would allow to determine if the field correlation 
length is much larger, much smaller, or compara- 
ble to IMpc. 

6. High energy Neutrinos 

6.1. GRB neutrinos, - 10 14 eV 

Protons accelerated in the fireball to high en- 
ergy lose energy through photo-meson interaction 
with fireball photons. The decay of charged pions 
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produced in this interaction, tt + — > fi + — > 
e + + v e + 77^ + fju, results in the production of 
high energy neutrinos [^3| . The neutrino spec- 
trum is determined by the observed gamma-ray 
spectrum, which is well described by a broken 
power-law, dN 1 jdE 1 oc E~^ with different val- 
ues of f3 at low and high energy 0]. The ob- 
served break energy (where (3 changes) is typi- 
cally E\ ~ lMeV, with f3 ~ 1 at energies below 
the break and /3 ~ 2 above the break. The inter- 
action of protons accelerated to a power-law dis- 
tribution, dN p /dE p oc E~ 2 , with GRB photons 
results in a broken power law neutrino spectrum, 
dN v jdE v oc E-P with /3 = 1 for E v < E h v , and 
j3 — 2 ioY E v > E h v . The neutrino break energy 
E h v is fixed by the threshold energy of protons for 
photo-production in interaction with the domi- 
nant ~ 1 MeV photons in the GRB Q, 

E b v w 5 x 10 u r2 00 (^/lMeV)- 1 eV. (5) 

The normalization of the flux is determined by 
the efficiency of pion production. As shown in 
Q , the fraction of energy lost to pion production 
by protons producing the neutrino flux above the 
break, E h v , is essentially independent of energy 
and is given by 

U U - 2 (^/lMeV)rf 00 Ai 10ms - 1 ' 

Thus, acceleration of protons to high energy in 
internal fireball shocks would lead to conversion 
of a significant fraction of proton energy to high 
energy neutrinos. 

If GRBs are the sources of UHECRS, then us- 
ing Eq. (ph and the UHECR generation rate given 
by Eq. (Q) with £ ~ 1, the expected GRB neu- 
trino flux is Q 

» 1.5 x 1CT 9 (|0 x 

min{l,£W^}-S^, (7) 
cm^ssr 

where v x stands for u^, and v e . 

The neutrino spectrum (Q) is modified at high 
energy, where neutrinos are produced by the 
decay of muons and pions whose life time ex- 
ceeds the characteristic time for energy loss due 
to adiabatic expansion and synchrotron emission 



|^,^5|,Q. The synchrotron loss time is deter- 
mined by the energy density of the magnetic field 
in the wind rest frame. For the characteristic pa- 
rameters of a GRB wind, synchrotron losses are 
the dominant effect, leading to strong suppression 
of v flux above ~ 10 16 eV. 

We note, that the results presented above were 
derived using the "A-approximation," i.e. assum- 
ing that photo-meson interactions are dominated 
by the contribution of the A-resonance. It has 
recently been shown , that for photon spectra 
harder than dN 1 /dE 1 oc E~ 2 , the contribution 
of non-resonant interactions may be important. 
Since in order to interact with the hard part of 
the photon spectrum, E 1 < E^, the proton en- 
ergy must exceed the energy at which neutrinos 
of energy E h v are produced, significant modifica- 
tion of the A-approximation results is expected 
only for E v ^ E h v , where the neutrino flux is 
strongly suppressed by synchrotron losses. 

6.2. Afterglow neutrinos, ~ 10 18 eV 

Protons are expected to be accelerated to > 
10 20 eV in both internal shocks due to variability 
of the underlying source, and in the reverse shocks 
driven into the fireball ejecta at the initial stage 
of interaction of the fireball with its surrounding 
gas, which occurs on time scale T ~ 10 s, compa- 
rable to the duration of the GRB itself. Optical- 
UV photons are radiated by electrons accelerated 
in shocks propagating backward into the ejecta, 
and may interact with accelerated protons. The 
interaction of these low energy, 10 eV-1 keV, pho- 
tons and high energy protons produces a burst of 
duration ~ T of ultra-high energy, 10 17 -10 19 eV, 
neutrinos [as indicated by Eq. (^)] via photo- 
meson interactions [£7J . 

Afterglows have been detected in several cases; 
reverse shock emission has only been identified 
for GRB 990123 |£l). Both the detections and 
the non-detections are consistent with shocks oc- 
curring with typical model parameters |p2| , sug- 
gesting that reverse shock emission may be com- 
mon. The predicted neutrino emission depends, 
however, upon parameters of the surrounding 
medium that can only be estimated once more ob- 
servations of the prompt optical afterglow emis- 
sion are available. 
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If the density of gas surrounding the fireball is 
typically n ~ 1cm -3 , a value typical to the inter- 
stellar medium and consistent with GRB 990123 
observations, then the expected neutrino inten- 
sity is fl7j 

2^F. ^ 10 10 



E„ 



-Y 



GeV 



\ 10 17 eV / cm 2 s sr ' 



(8) 



where = 1/2 for e° b - > 10 17 eV and (3 = 1 for 
e° b ' < 10 17 eV. Here too, v x stands for v^, and 
v e . The neutrino flux is expected to be strongly 
suppressed at energy > 10 19 eV, since protons 
are not expected to be accelerated to energy ^> 
10 20 eV. 

The neutrino flux due to interaction with re- 
verse shock photons may be significantly higher 
than that given in Eq. (0), if the density of 
gas surrounding the fireball is significantly higher 
than the value we have assumed, i.e. if n > 
Iran -3 . 

6.3. Implications 

The flux of ~ 10 14 eV neutrinos given in Eq. 
(0) implies that large area, ~ 1km 2 , high-energy 
neutrino telescopes, which are being constructed 
to detect cosmologically distant neutrino sources 
(see for review) , would observe several tens of 
events per year correlated with GRBs. The detec- 
tion rate of ultra-high energy, ~ 10 18 eV, after- 
glow neutrinos implied by Eq. (^) is much lower. 
The ~ 10 18 eV neutrino flux depends, however, 
on parameters of the surrounding medium which 
can be estimated only once more observations of 
reverse shock emission are available. 

One may look for neutrino events in angular 
coincidence, on degree scale, and temporal co- 
incidence, on time scale of seconds, with GRBs 
pTjf . Detection of neutrinos from GRBs could be 
used to test the simultaneity of neutrino and pho- 
ton arrival to an accuracy of ~ 1 s (~ 1 ms for 
short bursts), checking the assumption of special 
relativity that photons and neutrinos have the 
same limiting speed. These observations would 
also test the weak equivalence principle, accord- 
ing to which photons and neutrinos should suffer 
the same time delay as they pass through a grav- 
itational potential. With 1 s accuracy, a burst at 
100 Mpc would reveal a fractional difference in 



limiting speed of 10 -16 , and a fractional difference 
in gravitational time delay of order 10~ 6 (con- 
sidering the Galactic potential alone). Previous 
applications of these ideas to supernova 1987A 
(see |4j| for review), where simultaneity could 
be checked only to an accuracy of order several 
hours, yielded much weaker upper limits: of or- 
der 10~ 8 and 10 -2 for fractional differences in the 
limiting speed and time delay respectively. 

The model discussed above predicts the pro- 
duction of high energy muon and electron neu- 
trinos. However, if the atmospheric neutrino 
anomaly has the explanation it is usually given, 
oscillation to z/ T 's with mass ~ 0.1 eV [[30|, then 
one should detect equal numbers of v^s and f T 's. 
Up-going t's, rather than /i's, would be a dis- 
tinctive signature of such oscillations. Since i/ T 's 
are not expected to be produced in the fireball, 
looking for r's would be an "appearance experi- 
ment." To allow flavor change, the difference in 
squared neutrino masses, Ato 2 , should exceed a 
minimum value proportional to the ratio of source 
distance and neutrino energy A burst at 

100 Mpc producing 10 14 eV neutrinos can test for 
Am 2 > 10~ 16 eV 2 , 5 orders of magnitude more 
sensitive than solar neutrinos. 
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